Soft triaxial roto-vibrational motion in the vicinity of 7 = 7r/6 
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A solution of the Bohr collective hamiltonian for the /3— soft, 7— soft triaxial rotor with 7 ~ n/6 is 
presented making use of a harmonic potential in 7 and Coulomb-like and Kratzer-like potentials in f}. 
It is shown that, while the 7— angular part in the present case gives rise to a straightforward extension 
of the rigid triaxial rotor energy in which an additive harmonic term appears, the inclusion of the 
P part results instead in a non-trivial expression for the spectrum. The negative anharmonicities 
of the energy levels with respect to a simple rigid model are in qualitative agreement with general 
trends in the experimental data. 

PACS numbers: 21.60.Ev, 21.10.Ro 



The seareli for new solutions of the Bohr collec- 
tive hamiltonian has recently experienced a resurgence, 
mainly due to influential studies upon shape phase tran- 
sitions [1, 2]. The occurrence of analytically solvable 
models has served to define benchmarks for the anal- 
ysis of experimental nuclear spectra. As long as the 
quadrupole degree of freedom is considered, the nuclear 
Bohr-Mottelson collective model has a few well-known 
complete solutions that correspond to spherical or axially 
symmetric equilibrium intrinsic shapes. In some case, 
when the full solution is a demanding task, one has usu- 
ally resorted to the so-called rigid models in which one or 
both of the deformation variables, /3 and 7, is constrained 
to take a fixed value. The opposite situation, in which 
either the /? or 7 variable (or both) is not forced to take 
a fixed value, is called soft. 

The purpose of the present paper is to present a 
novel solution of the Bohr collective hamiltonian for 
the /?— soft, 7— soft triaxial rotor with 7 around 7r/6. 
Our model include /3 and 7 vibrations as well as rota- 
tions about a non-axially symmetric equilibrium intrinsic 
shape. 

Years ago Davydov [3] proposed a model for a rigid tri- 
axial rotor, described diagonalizing the hamiltonian for a 
rigid rotor with three different moments of inertia. The 
corresponding soft cases, that are described by the full 
Bohr hamiltonian with some V{P,'y) potential (that dis- 
plays a pocket whose minimum lies in /3o 7^ 0,70 ^ 0), 
has been only partially or numerically treated (see also 
[4] ) . Davydov discussed also a model with soft behaviour 
with respect to (3 [5] and a model with soft behaviour with 
respect to 7 [6]. These approaches give rather compli- 
cated expressions for the spectrum. The latter, in which 
the effective potential in 7 is expanded in a series and 
truncated neglecting terms of order higher than the sec- 
ond, has some similarity with the present approach. Our 
model, although restricted to 7 ^^ 7r/6, is easier to han- 
dle and gives in a straightforward fashion the full energy 
spectrum, while the cited Davydov's treatments give only 
the energies of a limited number of states. 



This study is partly born out from the necessity to 
overcome the reservations on the above models (see for 
example [7] ) about the neglected role of the fiuctuations 
of the dynamic moments of inertia and partly to show 
that the naive idea that the spectrum of a soft triaxial 
rotor will be the sum of a rigid rotor part plus a harmonic 
/3— and harmonic 7— vibrations is essentially incorrect. 
In the present case the solution of the 7— angular part 
is indeed a straightforward extension of the rigid triaxial 
rotor energy in which an additive harmonic term appears, 
but the inclusion of the /3 part gives rise to a peculiar 
combination of the various terms. 

A closely connected solution of the Bohr hamiltonian 
has recently appeared as a preprint [8]. The common 
point between this work and ours is the presence of a 
minimum in 7 = 30°. In the cited work a critical point 
symmetry, called Z(5), is introduced for the prolate to 
oblate shape phase transition. The authors discuss the 
solution and compare their predictions with experimental 
data in the Pt region, obtaining a good agreement. The 
main difference between the two approaches lies essen- 
tially in the expressions of the potentials. Their poten- 
tial may be written as mi(/3) -I- ^2(7), where a harmonic 
oscillator is used for the 7 variable, and where an infinite 
square well is used for the f3 potential, similarly to what 
has been done in Refs. [1, 2] for the so-called E(5) and 
X(5) symmetries, thus achieving an approximate separa- 
tion of variables. Anticipating the following discussion 
we mention that in the present solution a potential of 
the form (5) makes exact the separation of variables and 
Coulomb and Kratzer potentials are used in the f3 vari- 
able. 

Two very recent papers [9] deal with a similar subject: 
beside a detailed study on phase transition in nuclei in 
the framework of the interacting boson model, a sum- 
mary of different solutions of the Bohr hamiltonian is 
done and an interesting solution for a triaxial case with 
7 ^ 7r/6 is discussed. The author assumes a potential of 
the form F(/?,7) = m(/3) + ^D(3^ cos3j and follows an 
approach similar to the one used in rcf . [2] . 

The Schrodingcr equation with the Bohr hamiltonian 
may be written as 
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i/B*(/3,7,0,) - S*(/3,7,0,) 



(1) 



where i/s = r/3+T^+T™t + T/(/3, 7) [10] and the kinetic 
terms have the following expressions 
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where the Qk, are the projections of the angular momen- 
tum (L) in the intrinsic frame. The potential term is in 
general a function of /3 and 7. The conditions, that must 
be imposed on this potential term in order to separate 
the above differential equation, have been discussed ex- 
tensively in the literature (see for example [5, 11, 12]). 
For our purposes it will be sufficient to notice that when- 
ever the potential is chosen as 



F(/3,7) = Vm + ^ 



(5) 



the Schrodinger equation is separable [12] into the fol- 
lowing set of second order differential equations: 

(T,3 + Vm-E + ^)f{(3) = (6) 

(^/3\T^ + Trot)+V2{-/)-n^<^{l,9,) =0 (7) 

with *(/3,7, 6*,;) = /(/3)$(7, 6li). Notice that the second 
equation does not depend on /3. This separation of vari- 
ables invokes a constant, fi. Here we solve the problem 
of a /3— soft, 7— soft triaxial rotor with a potential that 
has a minimum located in /3 = /3o and 7 = 70 = 7r/6. 
The simplest approach to this problem is to choose a dis- 
placed harmonic dependence for the 7 variable in order 
to have a pocket around the minimum. The potential 
in the asymmetry variable must be an even function of 
7 because of the indistinguishability with respect to the 
change of sign of 7, and must also be a periodic func- 
tion to assure that the operation of relabeling of the axes 
does not change the shape of the ellipsoid, but only its 
orientation. This statement is obviously not fulfilled by 
our choice that is to be considered therefore only locally 
in the region around the minimum. We notice that a 
discrete rotation of mr/3 in the polar coordinate system 
(/?, 7) does not affect our conclusions. 

In the 7— rigid case with a situation of maximum asym- 
metry ( for 7o = 7r/6) two of the moments of inertia 
are 'accidentally' equal [13, 14], while the axes of the el- 
lipsoid have not equal length. The hamiltonian is thus 
axially symmetric around the intrinsic 1-axis, that may 
be taken as a quantization axis, as pointed out in ref. 
[4] . We will restrict to a small region around tt/6 to take 
advantage of this simplification. Before discussing how 
to solve the 7— soft case, we multiply the two equations 
above by 2Bm/h^ and we define reduced energy and po- 
tentials as e = 2B„iE/h? and Ui{z) = 2BjnVi{z)/h? 



with (i — 1,2). Lower case symbols denote reduced quan- 
tities and operators. The new set of 'reduced' equations, 
given the definitions U = {2Bm/1^)0^ Ti {i ~ 7, rot) and 
u) = 2B„n/h^, is 



-^|^'|+-(^^)-^V)^(^) = « («) 



(t^+trot+U2{l)-UJ^^{l,e,) 



^ 



(9) 



The spectrum is determined by the solution of the first 
differential equation in which uj, that it is found from the 
solution of the second differential equation, plays the role 
of the coefficient of a 'centrifugal' term and, as we will 
show, yields a non-trivial expression for the energy levels. 
Around 7r/6 we may set 7 = 7r/6 + x, with x a small 
expansion parameter. It is easily seen that in this case 
the rotational part of the Bohr hamiltonian, through use 
of trigonometric relations, becomes 
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where we have added and subtracted 4(5i. 

The first step to solve the (7, 6*^) part of the problem 
is to change variable in Eq. (9), using 7 = 7r/6 + x and 
introducing the following simplifications: 

sin 37 = cos 3a; ^ 1 cos a; ^ 1. 

We obtain, specifying the harmonic dependence of the 
potential, U2{x) — Cx^, and exploiting relation (10), 
the following equation 

i-^ + Cx^-oj+Q --^Ql\<^{x,e,) = (11) 



ffl 



where the second order differential operator Hi has been 
defined. The wavefunction $(a;,0i) may be written, fol- 
lowing [4], as 



y i6J(UL,o) ^^"-^'^^'^^^^^^-^^^'^+^'^^^^^'-^^^'^^ 

(12) 
where the angular part is written in terms of Wigncr 
functions labeled by the projection of the total angular 
momentum on the 1-axis, i?, that is a good quantum 
number, while the functions 77,1 ^l,r{x) are eigenfunc- 
tions of the one dimensional Schrodinger equation for the 
harmonic oscillator, Hirjn ^l,r{x) = 0. The index n^ is 
the quantum number associated with the vibrations in 



the 7 degree of freedom. The action of the operators in 
eq. (11) on the respective wavefunctions gives 

^L,fl,n, =^/C(2n^ + l) + L(L + l)-^i?2 (13) 

The first term corresponds to a harmonic 7— vibration, 
while the second and third account for the rotational en- 
ergy of the 7 = 7r/6 rigid triaxial rotor, reproducing the 
Meyer-ter-vehn formula [4] . 

The complete spectrum may then be obtained solving 
eq. (8) with some appropriate potential, ui(/3), as, for 
instance, in Ref. [1, 2, 11, 12, 15-18]. The eigenvalue 
equation in the (3 variable may in general be written as 



x"(/3)+(£-"i(/3)-^)x(/3) 



= 



(14) 



having imposed /(/3) = P~'^x(j3). When the poten- 
tial ui{(3) is taken to be of a [16, 18] Coulomb-like form 
—A/ (3, or of a Kratzer-like form —A/ (3 + B/f3'^, with the 
substitutions z = 2/3-^6, e = — e, k = Aji^^fe) and 
/i^ = 9/4 -I- B -I- w, we can recast equation (14) as the 
Whittaker's differential equation: 
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whose well-known solution is expressed in terms of Whit- 
taker's functions. As in [16, 18] the requirement of 
the correct asymptotic behaviour fixes the spectrum, 
introducing the n^ quantum number associated with 
/?— vibrations. The reduced eigenvalues are 



e{ny,nf3,L,R) = 



A^/A 



(^^9/4 + B+iJL,B.n-, + 1/2 + n/3] ' 



(16) 

where lo comes from eq. (13). Energies are usually rede- 
fined fixing the ground state to zero and using the energy 
of the first 2+ state as unit, namely 



e(n^,n/3,L,i?)-e(0, 0,0,0) 
e{n„np,L,R) - ,(0, Q, 2, 2) - e(0, 0, 0, 0) 



(17) 



In fig. (1) we displayed, as a function of B and for a 
fixed value of C = 1 , the lowest states of the ground state 
band {np = 0, n-^ = 0) and the lowest O"*" of the /3— band 
(np = 1, ri^ = 0). It may be seen that when _B — *■ 00 the 
typical rigid rotor energies are recovered from the limit 
of the ground state energy levels. The relative position of 
the /3— band (eventually tending to infinity), with respect 
to the ground state band, gives an idea of the degree of 
/3— softness. The negative anharmonicities evidenced in 
this figure (that are strictly valid for a triaxial nucleus 
with 7 ~ 30°) are qualitatively consistent with most ex- 
perimental data on triaxial nuclei where typically the 3"*" 
and 4+ states are found at an energy that is lower than 
the rigid rotor prediction. In fig. (2) we plotted the low- 
est 2+ state of the ground state band for reference and 




FIG. 1; Reduced energies of the lowest state of the /3— band 
(dashed line) and of a few lowest states of the ground state 
band (solid lines) as a function of B. The limits for the energy 
levels when B —> 00, that correspond to the rigid triaxial rotor 
energies, are reported in the right side. Here we fixed C = 1. 




FIG. 2: Reduced energies of the {J'^ = 2+,n-, = 1,71/3 = 0) 
state (solid line) as a function of B, for various values of the 
strength of the harmonic potential in 7, C. The first 2^ of 
the ground state band (dashed line) is reported for reference. 



the {J^ — 2'^,n^ = 1,71^ = 0) state for various values of 
the parameter C, that is the strength of the harmonic os- 
cillator in 7. A low value for C corresponds to a soft well, 
while a large value corresponds to a more rigid situation. 
When B ^ 00 the states of the 7— band tend to finite 
values, as far as C is finite. The Wilets-Jcan limit (com- 
plying with the choice of Coulomb and Kratzcr potentials 
in f3 [16]) is correctly obtained from formula (16) when 
the 7— potential becomes flat (C = 0). The reduced en- 
ergy spectra shown in both figures do not depend on the 
parameter A. This model has, from one side, a limited 
applicability because it is confined to some special tri- 



axial nuclei for which the minimum of the potential lies 
at 7 = 7r/6, but from the other side, it has a consider- 
able flexibility as far as the position of the excited bands 
is concerned. The interplay between the two parameters 
may be exploited to fit experimental data. Notice that B 
and C do not separately fix the position of the respective 
/3 and 7 bands, but they both take part in a non-trivial 
way to determine the energy levels. It is nevertheless 
very well known [19] that many predictions obtained for 
a 7— rigid model nearly agree with an axially symmet- 
ric 7— soft model whenever the "frms of the latter equals 
the rigid value of the former. This fact may increase the 
applicability of the model. 

The addition of 7 and /3 vibrations to the rigid triaxial 
rotor spectrum (around 7 = 30°) shows interesting fea- 
tures: the well-accepted naive idea that the spectrum of 
the soft rotor is approximatively the sum of a rigid rotor 
term plus separate vibrational terms in f3 and in 7 is not 



necessarily correct. Indeed we have shown that, while 
the solution of the 7— angular part of the problem gives 
a straightforward extension of the rigid rotor formula in 
which a simple harmonic term for the 7 degree of freedom 
appears, the full spectrum is rather a more complicated 
function that essentially depends on the choice of the po- 
tential in p. Other choices are possible as, for instance, 
a displaced harmonic oscillator or a Davidson potential 
[15; 17] in (3, that lead to exact solutions. The solution 
presented here may well serve as a good starting point 
to analyze energy spectra in various mass regions where 
triaxiality is expected to play a relevant role. 
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